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ABSTRACT 



1 Based on observations collected at the European Southern Observatory, Chile, for proposals 076. B- 
0646(A) and 077B-0738(B). 
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We present near-infrared (H- and K-band) integral-field observations of the 
circumnuclear star formation rings in five nearby spiral galaxies. The data, ob- 
tained at the Very Large Telescope with the SINFONI spectrograph, are used to 
construct maps of various emission lines that reveal the individual star forming 
regions ("hot spots") delineating the rings. We derive the morphological param- 
eters of the rings, and construct velocity fields of the stars and the emission line 
gas. We propose a qualitative, but robust, diagnostic for relative hot spot ages 
based on the intensity ratios of the emission lines Br7, He I, and [Fell]. Appli- 
cation of this diagnostic to the data presented here provides tentative support 
for a scenario in which star formation in the rings is triggered predominantly at 
two well-defined regions close to, and downstream from, the intersection of dust 
lanes along the bar with the inner Lindblad resonance. 

Subject headings: galaxies: nuclei - - galaxies: ISM - - galaxies: kinematics 
and dynamics — galaxies: individual(NGC 613, NGC5248, IC1438, NGC1079, 
NGC1300) 



INTRODUCTION 



In many spiral galaxies of early- and intermediate Hubble type (Sa-Sc), active star 
formation is organized in a ring-like structure. These star formation rings offer a unique 
opportunity to study massive star formation in external galaxies: t hey produce up to 2/3 of 
the bolometric luminosity of their host galaxies (e.g. NGC7469; iGenzel et al.lll995l ). and 
often contain a large fraction of the entire star formation activity of the galaxy. 

The general picture of why molecular gas assembles in a ring is well understood as a 
natural consequence of a non-axisymmetric gra vitational potential, nea r ly always due to the 



presence of a stellar bar or oval distortion (e.g . ICombes fc Gerinlll985l ; iKnapen et al.lll995 



Heller fc Shlosmanlll996l : lButa fc Combeslll996l ). Because of its dissipative nature, molecular 
gas accumulates around the radii at which the stellar orbits experience dynamical resonances 
with the rotating bar potential. Depending on the pattern speed of the bar and the rotation 
curve of the galaxy, there can be one or multiple such resonances. The high gas densities 
found in the rings, combined with a variety of excitation mechanisms such as ultra-violet 
radiation from young stars and mechanical shocks due to energetic outflows from massive 
stars and/or an active galactic nucleus (AGN) help reveal the physical state of the interstellar 
matter (ISM), be it molecular, atomic, or ionized gas. 



Besides being fascinating laboratories in their own right, star formation rings are impor- 
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tant also for the secular evolution of disk galaxies (IKormendy fe Kennicuttl 12004 ) . This is 
particularly true for the innermost of the dynamical resonances which is called either the "in- 
ner Lindblad resonance" (ILR) or - in cases where a compact massive object leads to an addi- 
tiona l dynamical resonance - the "nuclear Lindblad resonance" (NLR; iFukuda. Wada. fe Habe 
19981 ). Either of these resonances can produce gas rings with radii of a few hundred pc or 
less, depending on the enclosed mass, and the rotation speed of the galaxy disk. On these 
spatial scales, a number of processes can cause the gas to lose angular momentum and to 
subsequently flow to wards the nucleus. Exampl es for such processes include torques due to 
the stellar potential (iGarcia-Burillo et al.ll2005l ). dynamical frictio n between giant molecu- 



lar cl ouds that form within the ring due to se lf-gravity of the gas (IFukuda. Habe. fe Wada 



20001 ) . the formation of spiral density waves (jEnglmaier fe Shlosmanll2000l ). or mechanical 
energy released by star formation in the ring via stellar winds and/or supernova explosions. 

Understanding the gas dynamics and star formation processes of nucleaij^ rings in disk 
galaxies is therefore crucial for developing models of gas accumulation at the very nucleus of 
a galaxy and the evolution of any compact mass ive object (CMO) which can exist in the form 
of a nuclear star cluster (NC. lBoker et al.ll2002l ) and/or a supermassive black hole (SMBH). 



Whil e some theoretical models exist on the gas behavior around a CMO (e.g. IFukuda et al. 



20001 ). observational data to constrain these models are rare, mostly because of the limited 
spatial resolution of mm-observations. Only recently has it become possible to study the 
molecular ga s flows within a few tens of pc fro m the nucleus in a small number of nearby 



galaxies (e.g. ISchinnerer et ajj|2003j , 12009 . 120071 ) 



In order to increase the number of well-studied nuclear rings, we have begun a project 
to study the near-infrared (NIR) properties of five such objects, using the SINFONI integral- 
field spectrograph on the Very Large Telescope (VLT). This paper discusses the morpholo- 
gies, star formation rates, and kinematic properties of the rings. In a subsequent paper, 
we will make use of the spectroscopic information contained in the SINFONI data to in- 
vestigate in more detail the star formation process, stellar populations, and gas excitation 
mechanism(s) in individual galaxies, both in the rings and the galaxy nuclei. 

This paper is structured as follows. In §2J we describe the sample selection, observational 
details, and data reduction techniques common to all galaxies. The continuum and emission 
line morphology of the rings as well as the velocity fields of stars and gas for the individual 
galaxies are presented in £|3j We discuss the results in the context of competing models for 
the propagation of star formation in the rings in and summarize our analysis in S|5l 



2 Throughout this paper, the term "nuclear ring" is used to imply that it is the innermost (star forming) 
gas ring that can be resolved with the resolution of our data, typically a few hundred pc in diameter. 
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OBSERVATIONS AND DATA REDUCTION 



2.1. Galaxy Sample 



The general properties of the five galaxies discussed in this paper are summ arized in 
Table [H The objects were selected from the imaging survey of iKnapen et al.l (120061 ) because 
they show evidence for small-scale circumnuclear star formation rings within the SINFONI 
field of view, i.e. the central 8". All five galaxies are classified as barred (SB) or weakly 
barred (SAB). 



In most cases, the ring is obvious in the Ha images of IKnapen et al.l (120061 ) which have 
a spatial resolution comparable to that of the SINFONI data discussed in this paper. The 
only excepti on is NGC 1079 w hich requires HST resolution to unambiguously identify the 
nuclear ring (IMaoz et al.lll996f ). 



2.2. SINFONI observations 

SINFONI is an adaptive optics assisted, cryogenic NIR integral-field spectrograph, com- 
missioned at the VLT. It provides imaging spectroscopy of a contiguous, two-dimensional 
field of 64 x 64 spatial pixels in the wavelength range from 1.1 — 2.45 jim at a resolving power 



of 2 000 — 6000. Details a bout the instrument design can be found in lEisenhauer et al.l (120031 ) 
and Bonnet et al. J2004h 



The bulk of the data described in this paper were obtained in service mode during vari- 
ous nights between October 2005 and January 2006. One object (NGC 5248) was observed a 
semester later, in May/ June 2006. The SINFONI instrument was used without the adaptive 
optics module, and with its lowest magnification, resulting in slice widths of 0.25", a pixel 
scale of 0.125"/pixel, and a field-of-view of 8" x 8". This field size is well- matched to the di- 
ameters of the circumnuclear starburst rings under study, hence no mosaicing was necessary. 
We used the SINFONI configuration for simultaneous H+K spectra which yields a spectral 
resolution of 2000. 

Because in all cases, the host galaxy is much more extended than the SINFONI field of 
view, we obtained identical exposures of the empty sky background by "nodding" the tele- 
scope to point a few arcminutes away from the galaxy nucleus. During pipeline processing, 
these "sky frames" are subtracted from the "object frames", so that, at least in theory, the 
final data cubes are free of sky emission and dark current. 



For each galaxy, we had requested an on-source integration time of 2.5 hr, divided into 
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five observation blocks (OBs). Each OB consists of five 300s long on-source exposures, paired 
with corresponding observations of the empty sky in an ABBA nodding sequence. However, 
due to scheduling constraints, not all OBs were actually executed. In some cases (e.g. 
NGC 1079 and NGC 1300), only one or two OBs were executed. Although this clearly reduced 
the signal-to-noise ratio of the resulting data cubes, the emission line maps nevertheless yield 
interesting insight into the structure of the rings. In Table |2] we list the observation dates, 
OB identifiers, and the resulting total integration times for all galaxies discussed in this 
paper. 



2.3. Data reduction 

The reduced science products delivered by ESO for these datasets could not be used 
directly for our analysis. We identified two main reasons for this: 

a) First, the default settings of the ESO pipeline uses the flag OBJNOD.scale_sky 
= true in the Sl_REC_OBJNOD recipe to allow for scaling of the sky background before 
subtraction. This causes unphysical (often negative) flux levels in the resulting data cubes, 
especially in regions where the sky flux level is significant compared to the science object. 
It is thus important that this option is switched off when running the SINFONI pipeline on 
this type of raw data. 

b) Second, in the raw, two-dimensional detector frames, the sky emission lines show 
a small, but noticable shift in their position. This misalignment between sky and object 
frames causes imperfect removal of the sky lines in the resulting spectra, and is probably 
due to flexure in the optical train. 

While the first issue could be easily solved, the proper removal of the sky residuals 
required a more elaborate data reduction process. We used the ESO pipeline for SINFONI 
(version 1.3) to perform the basic data reduction (bias subtraction, flatfielding, wavelength 
calibration and extraction) of the individual frame pairs in each OB. For each pair, we im- 
proved on the sky subtraction by means of an IDDj software package described in iDavies 



(120071 ) and kindly provided to us by R. Davies. The algorithm allows to separately scale dif- 
ferent transition groups of telluric OH lines before subtraction which minimizes the residuals 
in the final spectra. We aligned all OBs to be centered on the galaxy nucleus (i.e. the peak of 
the continuum emission), thus accounting for telescope dithers between OBs. Moreover, we 
shifted all individual channel maps within an OB to correct the atmospheric refraction which 



3 http:/ /www. ittvis.com/ 
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causes small shifts with wavelength of the continuum peak within an exposure. Finally, we 
median-combined the aligned and sky-subtracted OB cubes of a given galaxy. 

The resulting "final" data cube of each OB contains about 2100 channel maps of the (8" x 
8") field of view with a spatial sampling of 0.125"/pixel and a spectral resolution of R~ 2000. 
The OB cubes are wavelength-calibrated and sky-subtracted, but at this stage neither flux- 
calibrated nor corrected for the absorption features introduced by the Earth's atmosphere. 
In order to correct for telluric absorption, we divided eac h data cub e by an atmospheric 



transmission model generated with the ATRAN package (iLordl Il992l ). We created model 
atmospheres for the altitude of the Paranal observatory and for a number of zenith angles 
spanning the range of our observations (0 deg — 30 deg). 

For flux calibration of the cubes, we used observations of early-type (around B5V) 
standard stars which are part of the routine nightly SINFONI calibration. For each night, we 
extracted the observed stellar spectrum (in counts/s) in a 5" aperture and - after correcting 
for telluric absorption as described below - divided it into the "true" intrinsic spectrum of 
the respective spectral type. This, in principle, yields the correction "spectrum" that must 
be multiplied into the data cubes obtained during the same night in order to obtain the 
correct SEP of the galaxy. The "true" stellar spectrum was obtained from the library of 



Pickles! (119981 ) after scaling to the near-infrared magnitudes listed in the Hipparcos catalo, 
For some stars, the exact spectral type was not available in the Pickles library, in which case 
we used the nearest available type (the difference was never more than one subclass). 

Performing this flux calibration on a night-by-night basis produced a set of galaxy SEDs 
that showed a relatively large scatter of about 15% rms. This scatter is larger than that 
achieved by using a single calibration spectrum for all OBs, regardless of when they were 
taken. This demonstrates that the accuracy of the flux calibration is limited by the system- 
atic uncertainty of the "true" spectra of the standard stars (and possibly by the accuracy 
of the telluric correction). We therefore decided to apply the same average calibration spec- 
trum to all individual OBs. This average calibration spectrum was obtained by fitting a 
fourth-order polynomial to the average inferred correction spectra of all standard stars. 

The systematic uncertainty of this approach should not be larger than the night-to-night 
variations of the flux calibration. We therefore estimate that our flux calibration is accurate 
to 15%. As a consistency check, we compared the continuum flux levels of the calibrated 
spectra in H- and K-band to those derived by the 2MASS survey. The agreement was better 
than 10%, thus confirming our error estimate. 



using the zero points of the UKIRT photometric system of 1.12 • 1(T 6 (4.07 • 1CT 7 ) erg cm 2 s 1 /im 1 for 
H-band (K-band). 
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As a last step in the data reduction, and after the telluric correction and flux calibration 
was applied to the individual OBs, they were merged to produce a fully reduced datacube 
for each galaxy. 



3. Near-Infrared Diagnostics: Stars, Gas, and Kinematics 

In this section, we present - for each of the five galaxies - a number of maps extracted 
from the SINFONI data cubes following the data reduction process just described. These 
maps contain a large amount of morphological and kinematic information on the stars and 
various components of the ISM (molecular and ionized gas, dust, etc.). The various emission 
lines, in particular, trace different components of the inter-stellar medium (ISM) and different 
excitation mechanisms so that even a qualitative comparison provides insight into the current 
physical conditions inside the starburst rings. 

Before we discuss each galaxy in more detail, we briefly summarize how each of the 
maps presented in Figures CD to [5] (from left to right and top to bottom) was constructed, 
what they signify, and how they can be interpreted: 



A map of the K-band continuum emission. In the NIR wavelength range, the con- 
tinuum emission is predominantly produced in red giants and supergiants, and the 
mass-to-light ratio of evolved stellar populations is nearly independent of stellar type. 
Therefore, the K-band morphology is an accurate tracer of the gravitational potential. 
Our K-band maps maps are constructed by simply adding up all channel maps in the 
wavelength range 2.0 — 2.4 jum, excluding those expected to contain strong line emis- 
sion, e.g. from Br7 or H 2 . The result should be a "clean" map of the stellar mass 
distribution. 



2. A 5 x 5 arcmin large-scale i?-band image of the galaxy from the Digital Sky Survey. 



A map of Ha emission in the central 8" x 8" taken from iKnapen et al.l (120061 ) . We 
include these maps here to allow a visual comparison to the SINFONI line maps, in 
particular Br7. However, the co mparison shoul d be m ade keeping in mind some compli- 
cations. First, the Ha images in lKnapen et al.l (120061 ) were constructed by subtracting 
a scaled I-b and image, used as a substitute for narrow-band continuum images. As 
discussed in iBoker et al.l (119991 ) . this method can cause imperfect continuum subtrac- 
tion in regions with large color deviations, in particular the galaxy nucleus. As can be 
seen by comparison to the SINFONI Br7 maps (which have a more robust continuum 
subtraction, see below), this is indeed the case for most of our sample. Secondly, Ha is 
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potentially more affected by dust extinction than Br7. Finally, the spatial resolution 
of the Ha maps is somewhat lower than that of the SINFONI data. Nevertheless, these 
Ha maps provide a valuable consistency check of our NIR line maps. Here, and in each 
of the following maps, we overplotted the K-band continuum contours from panel 1 for 
comparison. 

A [Fell] line map. The strongest emission line of [Fell] in the H- and K-bands is 
the a 4 Fg/2 — a A Dj/2 transition at 1.644 /iin. The most likely excitation mechanism is 
collisions with free electrons in hot gas with temperatures in the range 10 3 — 10 4 K. 
However, the ionization potential of [Fell] is only 16.2 eV. This means that in re- 
gions of fully ionized hydrogen, most of the Fe atoms are in higher ionization states. 
Therefore, the [Fell] emission is n ot expected to be strong in normal HII regions 



(IMouri. Kawara. fe Taniguchill2000l ). Instead, it will be strong in regions that i) have 
a high abundance of Fe in the gas phase, and ii) are only partially ionized. The two 
most likely mechanisms to create extended partially ionized zones are power-law pho- 
toionizatio n by X-rays (e.g. from an AGN) and fast shocks produced in supernova 



remnants (ILabrie fe Pritchet II2006I ). 



5. A He I line map, more specifically a map of the (n = 2) 1 P — 1 S transition at 2.058 fim 
which is the strongest He I line in the H- and K-bands. The ionization energy of 
atomic helium is 24.6 eV and thus markedly higher than that of hydrogen (13.6 eV). 
Line emission from He I is therefore expected to arise predominantly in the vicinity 
of the most massive (and hence youngest) stars. Interpretation of the He I emission 
requires detailed photo-ionization models, and quantitative analysis of t he temperature 



of th e hottest stars in an H II region is subject to large uncertainties (IDoherty et al. 



19951 ; iLumsden et al.ll200ll ). At least qualitatively, however, comparing the He I line 
strength to that of Br7 provides a handle on the relative ages of the youngest stellar 
clusters because the hottest stars will vanish fastest. We will use this simple diagnostic 
in the discussion of age gradients along the rings in §4.21 

A map of the Br7 line at 2.16 fim. The Br7 recombination line of the hydrogen atom 
is predominantly produced by photo-ionization in the Stromgren-spheres around O- or 
B-type stars. Although fast shock-fronts due to supernovae and/or strong stellar winds 
can also ionize hydrogen atoms, these mechanisms are related to the same population 
of OB stars, and thus they occur in the same locations^. A third potential ionization 



5 It is, in principle, possible to distinguish between ultra-violet photons and fast shocks as the cause of 
hydrogen ionization: the former produces low HII densities with a large filling factor (i.e. high column 
densities), while the latter produces high densities with a low filling factor. 
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mechanism is strong non-thermal emission from an active galactic nucleus (AGN). 
Outside the galaxy nucleus, however, we consider Br7 as a primary tracer of young 
stellar populations. 

The Br7 map, as well as the other emission line maps presented in Figures [T] to was 
constructed by summing up all channel maps containing line emission, and subtracting 
from each a continuum map created by averaging five channel maps on either side of the 
line. This is equivalent to the standard method of taking narrow-band images through 
filters that are centered on the line and the red and blue continuum, respectively. 

7. A H2 line map. The (1-0) SI transition of the H2 molecule at 2.122 fim is a tracer of 
warm molecular gas with temperatures around m 2000 K. There are three plausible 
mechanisms that can heat molecular gas to these temperatures: i) UV-fluorescence 
in photo-dissociation regions (PDRs), typically the surface l ayers of dense molecu- 
lar clouds in the vicinity of UV-bright star forming regions (e.g. lBlack fc yan Dishoeck 



1987 ). ii) collisional excitation of the H? molecules in fast shock fronts (e.g. 



19891 ) , and iii) X-ray illumination (IMaloney. Hollenbach fc Tielend 119961 ) . Using the 



Hollenbach fc McKee 



relative strengths of the numerous other rotational and vibrational transitions of the 
H 2 molecule present in the NIR spectral range, one can, in principle, distinguish be- 
tween these excitation modes. In the case of NGC 613, we will present this analysis 
in a forthcoming paper (IFalcon-Barroso et al.ll2008l . hereafter paper II). However, for 
the purpose of this paper, the H 2 (1 — 0) SI line maps merely serve to indicate the 
distribution of warm molecular gas. 

8. A map of the stellar velocity field, deriv ed using the penalized pixel- fitting (pPXF) 
method of ICappellari &: Emsellem I (120041) . W e made use of a subset of the library of 
theoretical spectra from lLancon et al.l (120071 ) adapted to our spectral resolution. This 
subset comprises stellar models with solar abundances for a range of effective tem- 
parature (T=2900-5900 K) and gravity (log(g)=0-2). In order to ensure reliable mea- 
surements of stellar kine matics, we spatially b i nned o ur datacubes using the Voronoi 
2D binning algorithm of ICappellari fc Copin I (120031 ). creating compact bins with a 
minimum signal-to-noise ratio (S/N) of ~ 30 per spectral resolution element. 

We fitted a non-negative linear combination of these models, convolved with a Gaus- 
sian line-of-sight velocity distribution (LOSVD), to the spectrum of each individual 
pixel over a wavelength range centered on the 2.29 /zm CO bands, masking any poten- 
tial emission lines. The fit also allows for a low-order Legendre polynomial in order to 
account for small differences in the continuum slope between the pixel spectra and the 
stellar templates. The best-fitting template mix was then determined by \ 2 minimiza- 
tion in pixel space. 
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9. A ma p of the Br7 velocity field, determined using the GANDALF code (jSarzi et al. 



20061 ) . This routine fits simultaneously both the stellar continuum and emission lines, 
treating the emission lines as additional Gaussian templates. In order to subtract the 
stellar continuum, we used the best template combination from the stellar kinematic 
fit described above. In general, the high signal-to- noise ratio of Br7 allows to reliably 
measure the gas kinematics in smaller spatial bins, which yields higher spatial resolu- 
tion than for the stellar kinematics, especially in the outer regions. However, the Br7 
kinematics can only be derived where the emission line is strong. In the figures, we 
only display velocities for regions with a S/N above 2, while regions below this S/N 
are left blank (gray areas). 

10. Finally, we present the H- and K-band spectrum of the galaxy nucleus, i.e. a circular 
aperture with 1" diameter centered on the K-band peak. Overplotted in red is the best 
stellar template fit from item 8. In the case of NGC613 (Fig. [1]), the emission lines as 
well as some prominent stellar absorption features are labeled for convenience. 



The quantitative analysis of the line fluxes over the entire SINFONI field as well as within 
an aperture encompassing the ring (as delineated in each of the Br7 maps) is summarized 
in Table HI In what follows, we comment in more detail on the individual galaxies. In 
particular, we describe those properties of the circumnuclear starburst rings that can be 
deduced from the images described above. 



3.1. NGC613 



NGC613 is a large, relatively nearby (d = 17.5 Mpc), strongly barred Sbc spiral. 
It is classified as a Seyfert galaxy in the NED database, although without a type desig- 
nation, and in fact, it is somewhat unclear on what basis this classification was made. 
The optical emission l i ne sp ectrum of the NGC613 nucleus is classified as "composite" by 
Veron-Cetty fc Veronl (119861 ). Furthermore, NGC613 emits neither detectable X-ray emis- 
sion in the 2 — 10 keV band , nor any high excitation MIR lines such as [OI V] or [Sill] that 
are usually found in AGN (jSatyapal. Sambruna. fc Dudiki 120041 ). In §3.71 we discuss that 
the nucleus is also quiescent in the NIR, in the sense that it does not ionize much of the 
surrounding gas. On the other hand, radio continuum observations show evid ence for an en 



erget i c outflow which has been i nterpreted as a jet emanating from an AGN (IHummel et al. 



19871 ; IHummel fc Jorsaterl Il992l ). The same radio observations also clearly reveal the star 
formation ring. 



The excellent signal-to-noise ratio of our SINFONI data for NGC613 allows a detailed 
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spectroscopic study of the star formation processes and stellar populations in the nucleus and 
nuclear ring of NGC 613 which we present in paper II. Here, we highlight a few interesting 
characteristics that are evident from the NIR maps alone. 



3.1.1. Continuum emission 



There is some dispute i n the literature about the shape of th e gr avitational potential in 
NGC 613. On the one hand. IJungwiert. Combes, k. Axon! (119971 ) and lEskridge et al.l (120021 ) 
have suggested t hat NGC 613 has a n uclear bar at the same position ang le (PA) as th e large- 
scale bar (110°; iHummel et al .11198 71 ) . However, based on HST imaging. lErwinl (120041 ) raised 
the possibility that the "nuclear bar" was in fact a spurious feature caused by the presence 
of an elliptical ring of star clusters. As can be seen from the following discussion, such a ring 
is indeed present, and the morphology of our K-band map in Figured] is clearly affected by 
the continuum emission from the ring. For this reason, and because of the limited SINFONI 
field of view we cannot decide whether the galaxy morphology underlying the ring shows a 
bar-shaped stellar mass distribution. 



3.1.2. Emission line morphology - evidence for an AGN outflow 



The Br7 emission shows a very clearly defined ring structure, composed of seven almost 
regularly spaced bright clumps that are the sites of current massive star formation. These 
"hot spots" are brightest along the southern half of the ring, while the northern half shows a 
well-defined "gap" at PA 30° which is also evid ent in the [Fe H i emiss ion map. This direction 
agrees well with tha t of th e radio jet found by IHummel et al.l (119871 ) and studied further by 
Hummel &: Jorsaterl (119921 ). It thus appears that a mechanical outflow from the AGN has 
disturbed the ring morphology. The outflow thus must intersect the plane of the disk. 



Additional evidence for such an outflow comes from the Ha image of iKnapen et al. 



( 120061 ) which shows two faint "streamers" that are closely aligned with the radio jet. More- 
over, a map of the velocity dispersion of the ionized gas (as measured by the line width of 
[Fell]) shown in Figure [6] clearly shows abnormally high dispersions along both sides of the 
outflow. The influence of the radio jet on the nuclear ring in NGC 613 will be addressed in 
more detail in paper II. 

The H 2 emission also shows the nuclear ring, but is strongest in the nucleus. The H 2 ring 
is not composed of distinct hot spots like the Br7 emission, but is more smoothly distributed. 
There are two emission peaks, found on opposite sides of the ring, at approximate PA of 90° 
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and 270°. If the H2 emission in the ring of NGC613 were purely caused by UV radiation, 
one might expect a spatial correspondence between H 2 and the ionized gas traced by Br7 
and/or He I. However, many of the Br7 hot spots are not bright in H 2 . Therefore, the H 2 
emission probably contains a non- negligible contribution from shock-heated molecular gas. 



3.2. NGC1079 



NGC 1079 is an isolated, grand-design spiral in the Eridanus supergroup (jBrough et al. 



20061 1 marked by a bright bulge and very faint outer disk. It represents a good example for 
a galaxy in which nearly the entire star f ormat ion activity takes place in the ring, as can be 
seen by the Ha image of iKnapen et al.l (120061 ) which shows an extremely quiescent galaxy 
disk 



The nuclear ring of NGC 1079 has first been identified by lMaoz et al.l (119961 ) from HST 
ultraviolet imaging and ground-based spectroscopy which confirmed that the gas rotation 
curve is consistent with a rotating ring. 



3.2.1. Continuum emission 



Within the rather small SINFONI field of view, the K-band map of NGC 1079 ap - 
pears smooth and featureless, in agreement with the H-band map of lJungwiert et al.l (119971 ) . 



Within the central 8", the K-band isophot es actually appear to have a somewhat smaller 
P. A. of ~ 75°, but this is also evident in the lJungwiert et al.l (119971 ) map. These authors find 
evidence for a secondary bar with P.A. 96° and major axis length of 17". It is possible that 
this bar governs the gas dynamics in the inner few hundred pc. Unfortunately, it is not pos- 
sible to determine the dust structure in NGC 1079 because no hi gh-resolution optic al/NIR 
imagery is available in the HST archive, except for the UV data of iMaoz et al.l (119961 ) which 
are however not suitable to derive extinction variations. 



3.2.2. Emission lines 

Despite its small angular size (radius ~ 1.5"), the star formation ring in NGC 1079 is 
easily resolved in our SINFONI maps. The line maps of Br7, H 2 , and [Fell] all show strong 
emission from the southern half of the ring, while the northern half is either quiescent or 
deeply obscured. The Br7 map shows two hot spots of roughly equal flux at P.A. f=a 150° and 
240°. Both are also seen in He I, but in this line the eastern hot spot is substantially stronger. 
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We detect little or no emission from H2 at the sensitivity limit of our data (note that the 
data set for NGC 1079 only consists of two OBs, i.e. only one hour on-source integration 
time). 



3.3. NGC 1300 



The SBbc spiral NGC 1300 is the prototype of a strongly barred disk galaxy, and 



has been extensively used as a testbed f or dynamical modelling (ILindblad fc Kristen 



Lindblad et al.lll997l : lAguerri et al. 



1996 



2000Jl The dust structure in the central kpc as revealed 



Martini et al. 



by the V-H color map of 
example of a grand design spiral. 



(j2003 ) does not resemble a ring, but rather is a good 



NGC 1300 shows no signs of nuclear activity, and lScarlatta et al.l (120041 ) find no evidence 
for a distinct stellar nucleus. Nevertheless, from a n HST/STIS study of the kinematics of 
the nuclear emission line gas, lAtkinson et al.l (120051 ) conclude that NGC 1300 hosts a SMBH 
of 6.6 ■ 10 6 M . 



3.3.1. Continuum emission 



As in the case of NGC 107 9, the continuum mo r pholo gy of NGC 1300 is very smooth, 
which has also been noted by IPerez- Ramirez et al.l (J2000|). There is no signature of the 



ring evident, indicating that in terms of stellar mass, the star formation ring is a small 
perturbation to the overall potential. 



3.3.2. Emission lines 

Of all galaxies in our sample, NGC 1300 is the one with the least integration time, 
and the SNR of the line maps, and in particular the velocity fields is therefore rather low. 
Nevertheless, the ring structure is clearly apparent, especially in Br7 and [Fell]. Both lines 
show a number of individual hot spots, with the two brightest Br7 clumps located at opposite 
ends of the ring along its major axis (P.A. « 125° and 305°). In contrast, the two strongest 
[Fell] clumps are located along the ring minor axis (P.A. fh 30° and 210°). The northern 
[Fell] clump is also bright in H 2 . 
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3.4. NGC5248 



The grand-design Sbc spiral NGC5248 is a "posterchild" ex ample for the impac t of a 
stellar bar on the gas dynamics and circumnuclear star formation. iJogee et al.l (j2002al ) have 
convincingly argued that the spiral structure is being driven by an extended, moderately 
strong stellar bar that has a deprojected ellipticity of 0.44 and a semimajor axis at, ar of 95" 



( 7. 1 kpc) . The in ner kiloparsec of NGC 5248 has been wel 



in em ission lines (Kennicutt et al.lll989l ; lLaine et al. 



2001 



studied with ground -based images 



Knapen et a 



1997 



Perez- Ramirez et al. 



tical ( Knapen et al. 2006 ) and NIR continuum (jElmegreen et al. 
20001 ) . These images have revealed a ~ 750 pc diameter star formation ring which, when 



200 6j) as well as op- 



seen at HST resolution, is clearl y resolved into numerous bright HII regions and young star 
clusters jMaoz et alJll996L l200lh . 



Near-infrared color maps have also revealed a nuclear dust spiral (ILaine et al.l 119991 ) 
that spans the radial range between 1"(75 pc) and 4"(300 pc). This structure likely forms 
the inward continuation of the grand-d esign spiral arm sy stem that is apparent in various 
tracers out to a radius of about 10 kpc (IJogee et al.ll2002al ). At its inside, the nuclear dust 
spiral appears to l e ad into a second H a-bright ring that has a radius of only 1.25" (100 pc) 
(ILaine et al.l l200ll ; iMaoz et al.l l200ll ) and thus forms one of the smallest such structures 
known. 

As with NGC 613, NGC 5248 is believed to harbor an AGN, and is classified as Sy2. 
However, at the sensitivity limits of the SINFONI data, we find no evidence for emission 
lines in the NIR spectrum of the NGC 5248 nucleus. 



3.4-1- Continuum emission 



The K-band m orphology in the center of NGC 5248 has been discussed in detail by 
Jogee et al.l (l2002al ). Our SINFONI map agrees well with the one presented in their Figure 8. 
In particular, both the "very weak oval feature of radius 3"" and the excess flux north of 



the nucleus are present in the two maps. As discussed by IJogee et al.l (j2002al ). the weak 
oval could be due to a disk-like component, a late-type bulge, or an unresolved nuclear bar. 
Our data cannot further constrain its origin, but we note that the stellar velocity field shows 
rather regular rotation over the SINFONI field, without any signs for kinematically distinct 
components. 
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3.4-2. Emission lines 



The two nested star formation rings of NGC5248 can be seen in the SINFONI Br7 
map. However, only the inner edge of the larger ring fall within the SINFONI field. The 
map clearly resolves the inner ring of NGC 5248. Th e morphology of th e inner ring agrees 
very well with the line maps obtained with HST by iMaoz et al.l (120011 ). in the sense that 
most Br7 hot spots can also be identified in Ha and/or Paa.The only region with significant 
He I flux is located in the inner ring at PA 170°. The strongest H2 emission is distributed 
over the inner 3", with a hot spot just south of the nucleus. We do not detect significant 
flux in the [Fe II] line within the SINFONI field. 



3.5. IC 1438 

At a distance of 33.8 Mpc, the poorly studied Sa spiral IC 1438 is the most distant galaxy 
in our sample, and with a radius of about 500 pc, the star formation ring is at the limit of 
what we consider "nuclear". 



3.5.1. Continuum emission 



Over the field of view of SINFONI, the K-band contours of IC 1438 are rather elongated 
outside the very nucleus. This i s unlikely to be a p rojection effect because IC 1438 has a 
rather low inclination of i = 32° (IKnapen et al.l 120061 ) . Indeed, the outer disk appears fairly 
round in the DSS image. On the other hand, the central 10 kpc of IC 1438 are domina ted 
by a strong stellar bar which is also evident in the I-band image of IKnapen et al.l (120061 ) . It 
is therefore likely that the gravitational potential in the center of IC 1438 is bar-shaped as 
well. It appears that the PA of the bar within the SINFONI field of view is higher than that 
of the large-scale bar. Whether this is the result of isophote twisting, or whether there is a 
"nested" secondary bar, requires additional data. 



3.5.2. Emission lines 

Due to the lower S/N of the SINFONI data, the Br7 map does not reveal the nuclear 
ring as clearly as the Ha map, mostly because the fainter, extended emission is not as 
obvious. Nevertheless, the two maps agree rather well, especially in the location of the "hot 
spots". Star formation is strongest at the south-east corner (PA 120°) where an extended 
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arc dominates the emission in both maps. This region also shows strong He I emission. Two 
smaller Br7 clumps at PA 0° and 260° are also apparent in the Ha map. There is also Br7 
emission in the nucleus, indicating star formation within the central ~ 50 pc. 

As in the other sample galaxies, the H2 emission is distributed much more smoothly, and 
thus outlines the ring structure more completely. There is some evidence for a connection 
between the ring and the nuclear region from an H 2 "ridge" at PA 260°. 



3.6. Kinematics of Stars and Gas 

In all five galaxies of our sample, the stellar velocity fields appear rather regular, despite 
the intense star formation within the nuclear rings. The PA of the kinematic major axes as 
listed in Table [I] was obtained by minimising the difference between the SINFONI velocit y 



fields and a bi(anti) symmetric model (for details, see Appendix C in lKrajnovic et al.ll2006l ). 
Within the uncertainties, the values are in good agreement with the PA of the photometric 
major axis, at least inside the SINFONI field of view. Unfortunately, the limited S/N of our 
data does not allow a reliable check for small misalignments and twists in the velocity fields 
that might be indicative for the presence of a stellar bar. 

Because Br7 is generally the most prominent emission line in our data, we used it for an 
estimate of the ionized gas kinematics. Our analysis of the Br7 velocity fields yields results 
similar to those derived from the stellar kinematics. In all cases, velocity gradients within 
the nuclear rings are smooth and the ionized gas has the same sense of rotation as the stars. 
Given the limited spatial coverage of the Br7 emission, we are unable to determine whether 
the ionized gas is kinematically misaligned with respect to either the stellar rotation field or 
the photometric major axis. 



In their analysis of NGC 1079, iMaoz et all (119961 ) noted that the gas inside the ring 



rotates faster than the stars at the same radius. Our observations confirm this result for all 
our objects: despite the limited spatial coverage of our Br7 velocity field, it is clear from 
the relevant panels in Figures [1] to [5] that the gas generally reaches higher velocities at both 
ends of the ring. This can be understood from the dissipative nature of the gas: because the 
kinetic energy of the gas is conserved, the rotation velocity increases while the dispersion 
decreases. As a consequence, any random motion in the z-direction is quickly dampened, 
and the gas settles in a thin disk. 

As mentioned in § 13.1.21 the case of NGC 613 warrants special attention, because the 
ionized gas is clearly disturbed by the radio jet. This is illustrated by the "twisting" of 
the isovelocity contours at the ends of the jet, and the increased velocity dispersion of the 
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surrounding gas as seen most clearly in the dispersion map of the [Fell] line presented in 
Figure El 



Two of our targets (NGC 613 and NGC 1300) are also present in the sample of lBatcheldor et al 



(120051 ) who derived stellar velocity fields for 23 spirals from integral-field observations over 
a comparable field of view, albeit at optical wavelengths. While their results for NGC 613 
agree well with ours, both in shape and amplitude of the stellar kinematics, there are signifi- 



cant d ifferences for NGC 1300. More specifically, the stellar rotation field of lBatcheldor et al. 



( 120051 ) (their Figure 3) appears "flipped" by 180° with respect to ours, i.e. it rotates in the 
opposite sense. We also note that their optical continuum map seems rotated by 90°. We 
believe that the maps presented in this study have the correct orientation, based on the 
comparison with the Ha map of iKnapen et all (120061 ) which agrees very well with our Br7 
map. 



3.7. Near-Infrared spectra of the nuclei 

In general, the nuclei (i.e. the central 1") of our sample galaxies appear rather quiescent. 
Except for NGC 613, which we will discuss in more detail below, the NIR spectra of the nuclei 
in Figures [1] to [5] show little, if any, evidence for NIR emission lines. There is faint nuclear 
[Fell] and H 2 emission apparent in the line maps of NGC 1300 (Figure [3]), as well as some 
Br7 emission in IC 1438 nucleus (Figure [5]). However, since these emission lines are not 
detected in the respective spectra, the apparent emission is likely a residual due to imperfect 
continuum subtraction. In Table HI we therefore list only upper limits for the emission line 
fluxes from the galaxy nuclei (except for [Fell] and H 2 in NGC 613). 

The lack of gaseous emission lines is not unexpected in a scenario in which gas accu- 
mulates at the nucleus over time until a critical density is reached. Star formation is then 
triggered, but can continue only until the gas supply is consumed or the energetic outflow 
from supernovae explosions disrupts and/or expels the gas, thus "quenching" star formation. 
This scenario is entirely consistent with observations, at least for late-type spirals: only about 
10% of nuclear star clusters currently show emissi on lines, although most of them harbo r 



a young stellar population less than lOOMyrs old (IWalcher et al.ll2006l ; iRossa et al.l 120061 ) . 
This can only be explained if star formation in galactic nuclei is episodic in nature, with a 
duty cycle of about 10%. 

This scenario appears to be supported by the case of NGC 613. The NIR spectrum 
of its nucleus shows several strong emission lines from the H 2 molecule, as well as strong 
emission from [Fell] at 1.644/xm. However, the Br7 recombination line is not detected at our 
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sensitivity limit, implying that star formation does not currently take place. This implies 
that while molecular gas is present in the nucleus, it has not sufficiently cooled to reach the 
critical density for star formation. The most likely excitation mechanism for both H 2 and 
[Fell] in the nucleus of NGC613 is the mechanical energy delivered by the AGN outflow. 
In the case of [Fell], this is different from the hot spots in the ring, where it is more likely 
produced by supernova remnants. In paper II, we will conduct a more detailed analysis of 
the physical state of the emission line gas in the nucleus of NGC613. 



4. Discussion: how does star formation proceed along a ring? 

The basic picture of star formation rings being the consequence of gas responding to 



a bar- shaped potential is generally accepted, and is nicely illustrated e.g. in lAthanassoula 



(119921 ): dynamical resonances associated with the non-axisymmetric rotating bar potential 
cause the orbiting gas to lose angular momentum and thus "spiral" towards the nucleus. In 
the process, the gas is concentrated in two elongated and curved dust lanes along the leading 
sides of the bar. In optical images, it appears as if the dust lanes connect the end of the bar to 
the nuclear ring. When observed at high angular resolution, the "ring" actually more closely 
resembles two tightly wound spiral arms. At the "contact points" between the dust lanes 
and the ring, the gas becomes less turbulent, and enters the almost circular orbits delineating 
the ring. While it is clear that there is abundant (molecular) gas throughout the ring, there 
is some debate about how and where star formation occurs. At least two mechanisms are 
plausible as illustrated in Figure [3, which we briefly describe in the following. 



4.1. "Popcorn" or "Pearls on a String"? 

In the first scenario, the gas accum ulates along the ring, until a critical density is 



reached. Presumably driven by turbulence (jElmegreenil 19941 ) . the ring then becomes unstable 



to gravitational collapse, and star formation is induced. In this model, either the entire ring 
structure begins to form stars at the same time, or individual hot spots collapse at random 
times and locations within the ring. Either way, there should be no systematic age sequence 
of hot spots along the ring. Because of its stochastic nature, we call this the "popcorn 
model" of star formation. 

In the second scenario, a short-lived, quasi-instantaneous "delta burst" of star forma- 
tion is induced in (and only in) a well-defined region along the ring where the gas density 
becomes sufficiently high to ignite star formation. The precise location of this "overdensity 
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region" (ODR) depends on the details of the potential, but it is often found close to, and 
somewhat downstream from where the gas enters the ring, as illustrated e.g. in Figure 15 of 
Regan fc Teubenl (l2003h . 



A young cluster that has formed at this location will continue its orbit along the ring, 
but star formation will cease as soon as the first supernova explosions expel the gas, i.e. 
after a few million years. Because the gas entering the ring is intrinsically clumpy, a series 
of short-lived starbursts will be triggered at the ODR. After star formation ends, the cluster 
will age passively as it follows its orbit. A series of such events will then produce a sequence 
of star clusters that enter the ring like "pearls on a string" . In this scenario, the star clusters 
should show a bi-polar age gradient around the ring, with the youngest clusters found close to 
the ODR, and increasingly older cluster ages in the direction of rotation, up to the opposite 
ODR. Because the age differences are simply due to the travel time along the ring, they 
should be consistent with the observed velocities of gas and stars in the ring, a question we 
briefly discuss in 14.31 

Evidence for such age gradients has been found in a few cas es. The most convi ncing; 
case is the nuclear ring in M 100, for which iRyder et all (120011 ) and lAllard et all (120061 ) used 
NIR and optical absorption indices and emission lines to determine the stellar ages of the 
hot spots in the star formation ring. ISmith et all (119991 ) have claimed evidence for sequential 
star form ation in the nuclear ring of NGC 7771 from the ratio of NIR to radio continuum. 
Moreover, iDavies et all (119981 ) have estimated the ages of individual hot spots in the ring of 
NGC 1068, based on the equivalent width of the Br7 emission. Although these authors do 
not specifically comment on this, their Table 2 shows hints for a bi-polar age gradient also 
in this star formation r i ng, a lthough the age differences appear too small to be significant. 
Finally, iMazzuca et all (120071 ) find evidence for such age gradients in 10 star formation rings 
out of a sample of 21. 



4.2. A new method for estimating hot spot ages 



In general, measuring the stellar cluster ages along the ring is difficult, even in a relative 
sense. The studies mentioned above have relied either on stellar absorption indices or equiv- 
alent width (EW) measurements of hydrogen recombination lines. The advantage of both 
diagnostics is that they are insensitive to extinction (but not to continuum emission from 
hot dust). According to the STARBURST99 models (ILeitherer et al.lll999l ). the equivalent 
width of hydrogen recombination lines produced in an instantaneous burst is more or less 
constant until about 3Myrs after the burst, and then drops steeply with time. It falls below 
the detection limit (about lA) of typical observations after only f^lOMyrs. Using e.g. the 
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EW of Ha is therefore limited to clusters in the age range 3 — 10 Myr, and any inferred age 
differences are small (few Myrs) by definition. 

Moreover, there is a fundamental difficulty with this kind of analysis: it requires accurate 
knowledge of the continuum emission from the hot spot itself, i.e. after subtraction of the 
underlying bulge and/or disk. In the centers of spiral galaxies, this separation between hot 
spot and bulge/disk is extremely difficult. The continuum emission from the bulge/disk is 
at least as bright as that from the hot spot itself and often much brighter, especially at 
NIR wavelengths where the light from evolved stars dominates. Because of the small field of 
view of integral field units, it is normally not possible from the data itself to disentangle the 
various components. Accurate continuum subtraction therefore often has to rely on uncertain 
assumptions on the galaxy structure. While the intrinsic symmetry of the bulge/disk will 
reduce the error for relative age measurements, estimating absolute ages of young clusters in 
circumnuclear star formation rings using absorption indices or emission line EWs is clearly 
fraught with uncertainty. 

Here, we use an alternative approach to obtain at least a qualitative sense for whether 
or not the young clusters in the ring show systematic relative age gradients. Instead of 
EWs, we use the flux of three emission lines that are prominent in the NIR spectra of star 
forming regions: He I, Br7, and [Fell]. Because line fluxes are independent of the underlying 
continuum level, they can be measured accuratehl without any knowledge of the disk/bulge 
structure. The assumption here is that the emission lines are produced only in the ring itself, 
i.e. that the underlying bulge/disk can be considered quiescent. Given that the line maps 
in Figs. [U to [5] show very little, if any, emission outside the ring proper, this seems a rather 
safe assumption for our targets. 

As mentioned in § El the Br7 and He I lines are both produced by photo-ionization in the 
Stromgren-spheres around O- or B-type stars. Because the He I line has a higher ionization 
energy than Br7, it requires the presence of hotter and more massive stars, and hence its flux 
falls off more rapidly after an instantaneous burst of star formation than that of Br7. The 
ratio of He I to Br7 flux can thus be used as a qualitative age indicator for cluster ages, in 
the sense that younger clusters still harbor hotter stars, and thus show more He I emission. 
We emphasize that that this is a qualitative statement only, and that we make no attempt 
to use the Hel/Br7 ratio to derive absolute ages, because the detai led interpretation of this 



line ratio is complicated (IDoherty et al.lll995t iLumsden et al.l 120011 ). 



6 None of the emission lines discussed here is significantly affected by possible stellar absorption features 
"underneath" the line. In fact, we did subtract the best-fitting mix of stellar templates before measuring the 
line fluxes listed in Table 21 However, this made little or no difference to obtaining fluxes from the spectra 
directly (i.e. without prior subtraction of the stellar continuum). 
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Assuming that all clusters have a similar initial mass function (IMF), and are sufficiently 
massive for statistical variations in stellar inventory to be insignificant, this diagnostic is 
useful for cluster ages up to about 10 M yr. Around that tim e, the supernova rate for an 
instantaneous burst reaches its maximum (ILeitherer et al.lll999l ). As described in §[31 the fast 
shock fronts associated with these supernovae are efficient in creating the extended partially 
ionized zones that produce strong [Fell] emission. Relative to He I and Br7 (which both 
have subsided significantly by this time), the [Fell] line becomes dominant for ages above 
« 10 Myr. 

The range of cluster ages prob ed by this approach is — 35 Myr, when the supernova 
rate drops almost instantaneously (ILeitherer et al.lll999l ). This time scale is significantly 
larger than that probed by EW measurements, and is well matched to the expected travel 
time of gas and star clusters around the ring: with a rotation velocity of 100 km s -1 , it would 
take about 15 Myr to complete a half orbit for a ring diameter of 1 kpc. As can be seen from 
Table [3j most of the rings in our sample are even smaller than this. The typical rotation 
velocities of the gas (which better probes the ring material) are of order 100 km s -1 (see the 
velocity maps of Br7 or the stars in Figures!]] to [5J and so it is entirely plausible that any age 
differences probed by these emission lines are indeed reflecting the travel times of clusters 
on their orbit along the ring. 

Finally, we note that for th ese relatively shor t orbit al timescales, the ring diameter can 
be regarded as constant. While iRegan &: Teubenl (120031 ) have shown that nuclear rings be- 
come smaller over time due to dissipation of angular momentum, this occurs over timescales 
of hundreds of Myrs, and thus can be ignored for the the present discussion. 



4.3. Application to the data 



In order to look for evidence of systematic age gradients as predicted by the "pearls on 
a string" scenario described in § 14.11 we first need to identify the location of any ODRs, if 
present. This is not a straightforward task, because maps and velocity fields of the molecular 
gas are difficult to obtain, especially for sources in t he southern hemisphere, and reach the 
required resolution only in a few nearby galaxies (e.g. ISchinnerer et al.ll2003l . 120061 ) . Lacking 
molecular gas observations, the best approach is to use hi gh spatial resolution maps of the 
dust morphology, as demonstrated by lMartini et al.l (120031 ). We have therefore searched the 
HST archive for optical images. Whenever such imagery was available, we identified dust 
lanes, their intersection with the ring, and any region(s) of high extinction along the ring 
that could plausibly be identified with an ODR. In cases where only one ODR was obvious, 
we have assumed that the second is located on the opposite side of the ring. 
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Figured] shows the case of NGC613. The left panel contains the HST/WFPC2 V-band 
(F606W) image of the central 15"(1.3kpc). Two prominent dust lanes are easily identified, 
and the inferred approximate locations of the ODRs are marked by star symbols. Note that 
the small size of the symbols certainly does not imply a positional accuracy, but merely 
indicates a ro ugh estimate. In fact, the ODRs can be quite extended as demonstrated in the 



simulations of iRegan fc Teubenl (120031 ). 



In both panels, the star formation ring is indicated by an ellipse with the ring parameters 
as listed in Table [3J Here, and for all other galaxies, the sense of rotation was determined by 
assuming that the outer spiral arms are trailing. In the right panel, we present a false-color 
(RGB) image constructed from the SINFONI emission line maps of He I, Br7, and [Fell], 
assigned to the blue, green, and red color channel, respectively. 

At least for the southern half of the ring, the color map appears to show the exact 
trend predicted by the "pearls on a string" scenario: starting at the ODRs, and moving in 
the direction of the gas flow, the color sequence is blue-green-red. As discussed in § I4.2[ 
this is exactly as expected if the clusters age passively along their orbits. The trend is less 
pronounced in the northern half of the ring, but here the picture is complicated by the strong 
mechanical interaction between the jet and the ring. It is possible that any cluster located 
within the "gap" at P.A. 30° has been stripped of its gas, and thus does not show up in the 
line maps, even though it might be relatively young. 

For a variety of reasons, the corresponding analyis for the other sample galaxies yields 
a somewhat less convincing picture. The case of IC 1438, presented in Figure [9j suffers 
from the fact that is has never been imaged with HST. In order to gain insight into the 
dust morphology, w e therefore used the B- and I-band ground-based images presented in 
Knap en et al. ( 2006 ) to create a B-I "pseudo" color map0. While the ring stands out as very 



blue in B-I (i.e. bright in Figure [H]), outside the ring the dust morphology is rather smooth. 
The only discernible structure is a faint region of higher extinction north of the nucleus 
which we interpret as the arm of a gas spiral. We also assume that one ODR is located close 
to the intersection of this spiral and the ring, and that the other lies opposite to it. With 
these - admittedly rather uncertain - assumptions, however, the picture for IC 1438 is similar 
to that for NGC613: He I peaks close to the eastern ODR, while [Fell] is strong in clusters 
that have traveled a good distance around the ring. 



7 The term "pseudo" is used here because no reliable magnitude scale can be attributed to the color map 



presen ted in Figure [H The reason is that a careful sky subtraction was not possible for the iKnapen et al 
120061) images because most galaxies completely filled the field of view of the CCD used. Nevertheless, the 
map is still valid in a relative sense, and can be used to identify "redder" regions that likely have higher 
extinction. 
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The case of NGC 1300 is presented in Figure [TUl Because only one OB was executed for 
this object, the quality of the SINFONI line maps is less than optimal. This is particularly 
true for the He I map which only shows tentative evidence for two emission peaks at P.A. 
135° and 315°. Intriguingly, however, this is exactly where one would place the ODRs based 
on the dust morphology derived from the HST/WFPC2 F606W map. Again, [Fell] emission 
peaks are located far from the assumed ODRs. 

Figure [TT] shows the false color map of the emission line morphology in the ring of 
NGC 1079. For this galaxy, there is neither HST imagery available, nor was it possible to 
identify any structure from a "pseudo" color map as described for IC 1438. We therefore are 
unable to estimate the location of the ODRs, and therefore cannot use this galaxy to test 
for a systematic age gradient along the ring. 

Finally, the case of NGC 5248 does not appear to be a valid test of the proposed scenario 
because it has two nested rings, the larger one of which falls mostly outside the SINFONI 
field of view (but can be recognized in the line maps of Figure H]). Because any gas falling 
onto the inner ring must first pass through the outer ring, and will be disturbed by star 
formation in the outer ring, the gas dynamics in NGC 5248 are likely very different from 
those in the other galaxies. This is possibly the reason why the dust morphology as shown 
in Figure [12] also appears different from that in the cases dicussed above. NGC 5248 harbor s 



a "flocculent" nuclear dust spiral without any clearly defined dust lanes (lLaine et al.ll200ll ). 
Because of the lack of dominant dust features, we were unable to determine the location of 
any putative ODRs. 

In summary, three of the five galaxies in our sample show line ratios that appear con- 
sistent with the "pearls on a string" scenario of how star formation progresses in nuclear 
rings. The data for the remaining two galaxies are incomplete, but not inconsistent with 
either model. 



4.4. Comparison to other studies 



Our findings lend some support to a picture in which young stars in nuclear starburst 
rings predominantly form in two distinct regions close to the "contact points" between the 
dust lanes and the ring. This is consistent with observations of the molecular g as distribution 



in the nuclei of barred galaxies which often shows a "twin peaks" geometry (e.g. lKenney et al. 
19921 ). This implies that cold molecular gas (i.e. the star formation fuel) exists only in two 



locations along the ring. Other studies claiming evidence for an a ge sequence of hot spots 
in nuclear rings have been mentioned in § 14.11 On the other hand, iMaoz et al.l (119961 ) have 
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found that the starburst rings in NGC 1079 and NGC 5248 experienced star formation over 
several 10 8 yr. This apparent contradiction might be explained by taking into account the 
dynamical timescales. As discussed in section 14.14 it takes as little as 15 Myr for a stellar 
cluster in the ring to rotate around half the ring. Therefore, several generations of stellar 
clusters will mix and overlap after typically 30 — 50 Myr. These older clusters are still 
prominent in broad band images (but not in the emission line maps discussed here). Using 
color diagnostics therefore might lead to a picture of "random" star formation in the ring, 
even if the "pearls on a string" scenario is correct. This highlights the importance of using 
age diagnostics that are well-matched to the orbital time scales. 

5. CONCLUSIONS 

Based on high-resolution NIR integral-field observations of five nuclear star formation 
rings, we have presented their emission line morphologies and velocity structure both in gas 
and stars. We have introduced a new method to derive relative hot spot ages along the rings 
using the relative strengths of the He I, Br7, and [Fell] lines. We employ this method to 
investigate the plausibilty of two competing scenarios for the way star formation is induced 
in nuclear rings, namely i) the "popcorn" model in which hot spots appear stochastically 
around the ring, and ii) the "pearls on a string" scenario in which star formation is triggered 
predominantly at two overdensity regions on either side of the ring. Only the latter predicts 
a well-ordered age sequence of hot spots along either half of a ring. 

The data presented in this study provide tentative support for the "pearls on a string" 
scenario, in that three out of five sample galaxies show some evidence for an age gradient of 
hot spots along the ring, while the remaining two galaxies have incomplete information and 
thus are consistent with either model. 

Of course, it might well be that star formation proceeds differently in some rings than 
in others. However, given that nuclear rings appear to form via a common mechanism (i.e. 
the gas response to a bar-shaped potential), it is not unreasonable to expect that they also 
follow a common path for inducing star formation. The small number of objects described 
here is clearly insufficient to provide reliable statistics, and similar studies of larger galaxy 
samples are needed to decide whether a particular mechanism governs the star formation in 
the majority of nuclear rings. Be that as it may, the proposed method of using NIR line 
ratios to estimate relative ages of young star clusters has been demonstrated to be a powerful 
tool for studies along this line. 
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(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


Galaxy 


Type 


Activity 


D 


scale 


M B 


i 


PA 


PA/jj n 


PA b ar 


Ibar 








(Mpc) 


(pen 


(mags) 


(°) 


(°) 


(°) 


(°) 


(") 


NGC613 


SB(rs)bc 


Sy 


17.5 


84.8 


-20.53 


41 


120 


120 


127 1 


> 59 1 


NGC1079 


SAB(r'l)a 




16.9 


81.9 


-18.83 


52 


87 


65 


122 1 


32 1 


NGC1300 


SB(s)bc 




18.8 


91.1 


-20.42 


49 


106 


116 


100 2 


94 2 


NGC5248 


SB(rs)bc 


Sy2 HII 


22.7 


110.1 


-21.07 


44 


110 


105 


105 3 


95 3 


IC 1438 


SAB(r)a 




33.8 


163.9 


-20.08 


32 


145 


158 


n.a. 


n.a. 



Table 1: Global parame ters of the five sample galaxies, obtained from the RC3 
(jde Vaucouleurs et al.l Il99ll ) unless otherwise indicated. Tabulated are the galaxy name 
(col. 1); the morphological type (from NED; col. 2); nuclear activity class (NED; col. 3); 
distance D in Mpc from Nearby Galaxies Catalog (lTullyill988l . col. 4); i mage scale in parsec 
per arcsec, as derived from the distance (col. 5); absolute blue magnitude (jTulfylll988l . col. 6); 
inclination % as derived from the ratio of the major to the minor isophotal diameter (col. 7); 
position angle (PA) of the disk photometric major axis (col. 8); PA of the kinematic major 
axis of the stellar velocity field, derived from this work (see § 13.61 c °l- 9); PA and semi-major 
axis length of the main stellar bar ( cols. 10 fc 11). Referenc es fo r the bar paramete rs are 
as follows: 1) ljungwiert etaD Jl997h . 2) iKnapen et aD J2002h. 3) Ijogee et al.l d2002bh . For 
IC 1438, the inclination estimate in Col. 8 is from iGarcfa-Gomez et al.l (120021 ). while no 
values are available for the bar parameters. 
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Table 2. SINFONI observations 



Galaxy 


tint 


OB 


Obs. date 


Seeing a ["] 






zloloo 


17.11.2005 


n r 
(J.D 






216157 


13.11.2005 


0.5 


NGC 613 


125 


216158 


23.10.2005 


0.5 






216159 


13.11.2005 


0.5 






216180 


27.11.2005 


0.6 


NGC 1079 


50 


216161 


17.11.2005 


0.6 






216162 


23.10.2005 


0.5 


NGC 1300 


25 


216168 


24.1.2006 


0.6 






216184 


10.5.2006 


0.7 






216187 


15.5.2006 


0.6 


NGC 5248 


125 


216188 


31.5.2006 


0.7 






216189 


1.6.2006 


0.5 






216190 


5.7.2006 


0.5 






216147 


22.10.2005 


0.7 






216169 


20.10.2005 


0.6 






216170 


19.10.2005 


0.5 


IC 1438 


150 


216171 


19.10.2005 


0.5 






216176 


20.10.2005 


0.6 






216278 


22.10.2005 


0.7 



a The atmospheric seeing was estimated from the 
FWHM of the K-band point spread function, mea- 
sured from standard star observations during during 
the respective nights. These were normally scheduled 
within one hour from the galaxy exposures. 
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Table 3. Summary of ring morphologies 



Galaxy 


D [pc] 


PA[deg] 


e 


NGC613 


603 


113 


0.40 


NGC 1079 


238 


72 


0.29 


NGC1300 


666 


117 


0.13 


NGC 5248 


338 


102 


0.39 


IC 1438 


1074 


81 


0.15 
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Table 4. Summary of emission line fluxes. 



Galaxv 


He I 




H 9 


[Fe III 


NGC6f3 










- ring 


46.5 


136.5 


80.5 


184.2 


- nucleus 


<0.5 


<1.1 


6.1 


11.8 


NGC 1079 










- ring 


6.2 


14.5 


0.8 


3.9 


- nucleus 


<0.2 


1.5 a 


<0.2 


<0.2 


NGC 1300 










- ring 


2.2 


19.6 


22.5 


15.5 


- nucleus 


<0.2 


<0.2 


<1.1 


<1.1 


NGC 5248 










- ring 


1.2 


4.5 


8.4 


<0.2 


- nucleus 


<0.2 


<0.2 


<1.0 


<0.2 


IC 1438 










- ring 


10.2 


13.8 


5.2 


10.6 


- nucleus 


<0.2 


<1.0 


<0.2 


<0.2 


Note. 


All 


fluxes 


are in 


units of 



10 15 erg cm 2 s 1 . Typical uncertainties 
are about 20% . 

a the central 1" aperture is likely "contam- 
inated" by flux from the nuclear ring. 
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Fig. 1. — Near-infrared morphology of NGC613. Please see § [3] for a description of the 
individual maps. Black contours delineate the K-band continuum from the top-left panel. 
The white ellipses in the Br7 map mark the ring aperture used for the line flux measurements 
listed in Table HI In all panels, north is up and east is to the left. 




Fig. 2. — As Figure CD, but for NGC 1079. See Figured] for identification of spectral features. 




Fig. 3. — As Figure CD, but for NGC 1300. See Figured] for identification of spectral features. 




Fig. 4. — As Figured! but for NGC5248. See Figure[T]for identification of spectral features. 
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Fig. 5. — As Figured], but for IC 1438. See Figured] for identification of spectral features. 
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Fig. 6. — Dispersion map of t he [Fell] emission line in N GC613. Note the line broadening 
along the radio jet at PA 30° (IHummel fc Jorsaterlll992l ) 
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"PcpwrnT "Pwris gn a siring:" 

Fig. 7. — Illustration of two candidate scenarios for star formation in a nuclear ring. Dark 
grey areas denote dense, cold gas that is conducive to star formation. The various star 
symbols denote young stellar clusters, their colors signify the cluster age in the sense that 
age increases from blue to green to red. A clear age sequence is expected only in the "pearls 
on a string" scenario (see § 14.11 for discussion). 



Fig. 8.— Left: Archival HST/WFPC2 F606W map of the inner 15" x 15" of NGC613. An 
ellipse with the parameters listed in Table [3] is overplotted in order to outline the assumed 
flow path of matter in the star formation ring. Also indicated are the assumed locations 
of overdensity regions (ODRs, marked by star symbols) where infalling gas is compressed. 
The sense of rotation is indicated with an arrow. Right: false color map of the emission 
line morphology in the inner 8" x 8" of NGC613. Here, the blue channel represents He I, 
green Br7, and red [Fe II] emission. Note that - at least in the southern half of the ring - 
the hottest stars (as marked by He I emission in blue) are located immediately downstream 
from the contact points. In contrast, the supernovae (as traced by the [Fell] emission in 
red) occur predominantly further along the ring, with Br7 emission dominated the region 
in between. This suggests that star formation is triggered at the ODRs, and that the star 
forming clumps then age passively along their orbit. 



Fig. 9. — Left: relative (i.e. uncalibrated) B-I color map of the inner 20" x 20" of IC 1438. 
An ellipse with the ring parameters as listed in Table [3] is overplotted in order to outline 
the assumed flow path of matter in the star formation ring. Also indicated are the assumed 
locations of ODRs (star symbols) where infalling gas is compressed. The sense of rotation 
is indicated with an arrow. Right: false color map of the emission line morphology in the 
inner 8" x 8" of IC 1438 as described in Figure [HI 
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Fig. 10.— Left: Archival HST/WFPC2 F606W map of the inner 15" x 15" of NGC 1300. An 
ellipse with the parameters listed in Table [3] is overplotted in order to outline the assumed 
flow path of matter in the star formation ring. Also indicated are the assumed locations of 
ODRs (star symbols) where infalling gas is compressed. The sense of rotation is indicated 
with an arrow. Right: false color map of the emission line morphology in the inner 8" x 8" 
of NGC 1300 as described in Figure E 




Fig. 11. — False color map of the emission line morphology in the inner 8" x 8" of NGC 1079 
as described in Figure [HJ The sense of rotation is indicated with an arrow. 



Fig. 12.— Archival HST/WFPC2 F814W map of the inner 15" x 15" of NGC5248. An 
ellipse with the parameters listed in Table [3] is overplotted in order to outline the assumed 
flow path of matter in the inner star formation ring. The sense of rotation is indicated with 
an arrow. We were unable to even tentatively identify distinct ODRs in this object. 



